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Abstract 

We report the regions where a signal-to-noise ratio (SNR) gain exceeding unity ex- 
ists in a parallel uncoupled array of identical bistable systems, for both subthreshold 
and suprathreshold sinusoids buried in broadband Gaussian white input noise. Due to 
independent noise in each element of the parallel array, the SNR gain of the collective 
array response approaches its local maximum exhibiting a stochastic resonant behavior. 
Moreover, the local maximum SNR gain, at a non-zero optimal array noise intensity, 
increases as the array size rises. This leads to the conclusion of the global maximum 
SNR gain being obtained by an infinite array. We suggest that the performance of infi- 
nite arrays can be closely approached by an array of two bistable oscillators operating 
in different noisy conditions, which indicates a simple but effective realization of arrays 
for improving the SNR gain. For a given input SNR, the optimization of maximum 
SNR gains is touched upon in infinite arrays by tuning both array noise levels and an 
array parameter. The nonlinear collective phenomenon of SNR gain amplification in 
parallel uncoupled dynamical arrays, i.e. array stochastic resonance, together with the 
possibility of the SNR gain exceeding unity, represent a promising application in array 
signal processing. 

1 Introduction 

The past decade has seen a growing interest in the research of stochastic resonance (SR) phe- 
nomena in interdiscipHnary fields, involving physics, biology, neuroscience, and information 
processing. Conventional SR has usually been defined in terms of a metric such as the output 
signal-to-noise ratio (SNR) being a non-monotonic function of the background noise intensity, 
in a nonlinear (static or dynamic) system driven by a subthreshold periodic input |TI. For 
more general inputs, such as non-stationary, stochastic, and broadband signals, adequate SR 
quantifiers are information-theoretic measures ^E]. Furthermore, aperiodic SR represents 
a new form of SR dealing with aperiodic inputs O . The coupled array of dynamic elements 
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Figure 1: A parallel array of N archetypal over-damped bistable oscillators. Each oscillator 
is subject to the same noisy signal but independent array noise. In this paper, we call ^(t) 
the input noise and r]i(t) the array noise. 



0011117] and spatially extended systems jSj have been investigated not only for optimal noise 
intensity but also for optimal coupling strength, leading to the global nonlinear effect of spa- 
tiotemporal SR By contrast, the parallel uncoupled array of nonlinear systems gives rise 
to the significant feature that the overall response of the system depends on both subthreshold 
and suprathreshold inputs PITUliniiniiniinillSllinillZlIIHlIini- In this way, a novel form 
of SR, termed suprathreshold SR ^21 1 attracted much attention in the area of noise-induced 
information transmissions, where the input signals are suprathreshold for the threshold of 
static systems or the potential barrier of dynamic systems ^2 El El El El El E] • In 
addition, for a single bistable system, residual SR (or aperiodic SR) effects are observed in 
the presence of slightly suprathreshold periodic (or aperiodic) inputs 201 ■ 

So far, the measure most frequently employed for conventional (periodic) SR is the SNR 
Pill Eg. The SNR gain defined as the ratio of the output SNR over the input SNR, also 
attracts much interest in exploring situations where it can exceed unity |28l 1291 1801 1811 1191 
E2ESII21ElESlEniE3EHlEH!- within the regime of validity of linear response theory, it 
has been repeatedly pointed out that the gain cannot exceed unity for a nonlinear system 
driven by a sinusoidal signal and Gaussian white noise UHl I24L 1251 1261 177| . However, beyond 
the regime where linear response theory applies, it has been demonstrated that the gain 
can indeed exceed unity in non-dynamical systems, such as a level-crossing detector |2H1; a 
static two-threshold nonlinearity E] , and parallel arrays of threshold comparators or 

sensors ^102 EH); and also in dynamical systems, for instance, a single bistable oscillator 
13 1^ ESI EOl El EHl EHIj a non-hysteretlc rf superconducting quantum interference device 
(SQUID) loop 0n|) and a global coupled network "5^. 

A pioneering study of a parallel uncoupled array of bistable oscillators has been performed 
with a general theory based on linear response theory |2H1, wherein the SNR gain is below 
unity. Recently, Casado et al reported that the SNR gain is larger than unity for a mean- 
field coupled set of noisy bistable subunits driven by subthreshold sinusoids 0. However, 
each bistable subunit is subject to a net sinusoidal signal without input noise. The conditions 
yielding a SNR gain exceeding unity have not been touched upon in a parallel uncoupled array 
of bistable oscillators, in the presence of either a subthreshold or suprathreshold sinusoid and 
Gaussian white noise. In practice, an initially given noisy input is often met, and a signal 
processor operating under this condition, with the feature of the SNR gain exceeding unity, 
will be of interest [221 ESI- The SNR gain has been studied earlier in the less stringent 
condition of narrowband noise |41| . In the present paper, we address the more stringent 
condition of broadband white noise and the SNR gain achievable by summing the array 
output, wherein extra array noise can be tuned to maximize the array SNR gain. As the 
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array size is equal to or larger than two, the array SNR gain follows a SR-type function of the 
array noise intensity. More interestingly, the regions where the array SNR gain can exceed 
unity for a moderate array size, are demonstrated numerically for both subthreshold and 
suprathreshold sinusoids. Since the array SNR gain is amplified as the array size increases 
from two to infinity, we can immediately conclude that an infinite parallel array of bistable 
oscillators has a global maximum array SNR gain for a fixed noisy sinusoid. For an infinite 
parallel array, a tractable approach is proposed using an array of two bistable oscillators, in 
view of the functional limit of the autocovariance function jl^l- We note that, for obtaining 
the maximum array SNR gain, the control of this new class of array SR effect focuses on 
the addition of array noise, rather than the input noise. This approach can also overcome a 
difficult case confronted by the conventional SR method of adding noise. When the initial 
input noise intensity is beyond the optimal point corresponding to the SR region of the 
nonlinear system, the addition of more noise will only worsen the performance of system 
|42| . Finally, the optimization of the array SNR gain in an infinite array is touched upon by 
tuning both an array parameter and array noise, and an optimal array parameter is expected 
to obtain the global maximum array SNR gain. These significant results indicate a series of 
promising applications in array signal processing in the context of array SR effects. 



2 The model and the array SNR gain 

The parallel uncoupled array of N archetypal over-damped bistable oscillators is considered 
as a model, as shown in Fig. ^ Each bistable oscillator is subject to the same signal- 
plus-noise mixture s(t) + S,{t), where s{t) = A sin(27rt/Ts) is a deterministic sinusoid with 
period Tg and amplitude A, and ^(t) is zero-mean Gaussian white noise, independent of s{t), 
with autocorrelation {£,{t)£,{0)) ~ D^S{t) and noise intensity D^. At the same time, zero- 
mean Gaussian white noise r]i(t), together with and independent of s{t) + £,{t), is applied to 
each element of the parallel array of size N. The N array noise terms f7i(i) are mutually 
independent and have autocorrelation {rji{t)rji{Q)) = DriS{t) with a same noise intensity 
. The internal state Xi (t) of each dynamic bistable oscillator is governed by 

for i = 1, 2, . . . , A^. Their outputs, as shown in Fig. ^ are averaged and the response of the 
array is given as 

1 ^ 

1=1 

Here, the real tunable array parameters Tq and Xh are in the dimensions of time and ampli- 
tude, respectively 21 . We now rescale the variables according to 

X^{t)/Xb X,{t), A/Xb ^ A, t/Ta ^ t, T,/Ta ^ T„ D^/{TaX^) ^ D^, D^/{TaX^) -> D^, 

(3) 

where each arrow points to a dimensionless variable. Equation is then recast in dimen- 
sionless form as, 

^ = X, (t) - xl it) + s{t)+ at) + it) . (4) 

Note that s(t) is subthreshold if the dimensionless amplitude A < Ac = 2/^/27 w 0.385, 
otherwise it is suprathreshold ■ 

In general, the summed output response of arrays y{t) = ^i(^) ^ random 

signal. However, since s(t) is periodic, y{t) will in general be a cyclostationary random signal 
with the same period Tg PJ- A generahzed theory has been proposed for calculating the 
output SNR According to the theory in the summing response of arrays y{t), at 
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any time t, can be expressed as the sum of its nonstationary mean E[y{t)] plus the statistical 
fluctuations y{t) around the mean E[y{t)], as 

y{t)^m+E[y{t)]. (5) 

The nonstationary mean E[y{t)] — {i/N)^^^^E[xi{t)] is a deterministic periodic function 
of time t with period Ts, having the order n Fourier coefficient 

Y,, = (^E[y{t)]exp{^t27T^)y (6) 

where {•■■) — (l/7s) J^' ■ ■ ■ dt. For fixed t and r, the expectation E[y{t)y{t + r)] is given by 

E[y{t)y{t + r)] = E[y{t)y{t + r)] + E[y{t)]E[y{t + r)]. (7) 

Then, the stationary autocorrelation function Ryyir) for y{t) can be calculated by averaging 
E[y{t)y{t + r)] over the period Tj, as 

Ryyir) - {E[y[t)y{t + T)]) 

= {E[y{t)y{t + r)]) + {E[y{t)]E[y{t + r)]) 

= Cyy(r) + (i?[y(t)]i?[y(t + T)]), (8) 

with the stationary autocovariance function Cyy(r) of y{t). The power spectral density Pyyiy) 
of y{t) is the Fourier transform of the autocorrelation function i?yy(T) 

/ + 00 
i?yy(T)exp(-«27rz/T)dT 

= ^[Cyy(r)]+ YnTjiv^^). (9) 

n— — oo ^ 

It is seen that the power spectral density Pyy{v) is formed by spectral lines with magnitude 
at coherent frequencies n/Ts, superposed to a broadband noise background repre- 
sented by the Fourier transform of Cyy(r). Note that E\y{t)y{t)\ = var[y(t)] represents 
the nonstationary variance of y{t)^ which, after time averaging over a period T^, leads to 
Cyy(O) = (var[y(t)]), the stationary variance of y{t). The deterministic function Cyy(r) can 
thus be expressed as 

Cyy(T) = (var[y(t)])Mr), (10) 

where the correlation coefficient /i(t) is a deterministic even function describing the nor- 
malized shape of Cyy(T), having a Fourier transform JF[/i(r)] = H{v). The power spectral 
density of Eq. Q can then be rewritten as 

PyyH = (var[y(t)])i/M+ i^«0(^ - f )• (H) 

n— — oo ^ 

The output SNR is defined as the ratio of the power contained in the output spectral line at 
the fundamental frequency 1 /Tg and the power contained in the noise background in a small 
frequency bin AB around l/Tg, i.e. 

^""^^^/^"^ = (var[y(t)])i/(l/T,)Ai?- ^^^^ 

In addition, the output noise is a Lorentz-like colored noise with the correlation time Tr 
defined by 

h{\T\ > Tr) < 0.05. (13) 
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Figure 2: Numerical output behaviors of a single bistable oscillator for A — 0.4 at four 
representative rms amplitudes a^. (a) The nonstationary mean E[y(t)]. (b) The stationary 
autocorrelation function Ryyir). (c) The stationary autocovariance functions Cyyir). Here, 
CTf = 0.6, 1.1, 1.8 and 3.0 correspond to curves a, f3, 7 and n. (d) The theoretical input 
SNR (solid line) of Eq. H14(l as a function of cr^. The numerical input SNR i?in (* but almost 
indistinguishable) and output SNR i?out (+) are also plotted. 



In the same way, the periodic sinusoidal input s{t) = yl.sin(27rt/rs) has total power /2 
and power spectral density + l/Tg) + 5{i^ — l/r5)]/4 in the context of bilateral power 

spectral density |31| . Here, the signal-plus-noise mixture of s(t) + ^(t) is initially given, and 
the theoretical expression of input SNR can be computed as 



_ AV4 
D^AB ~ ajAtAB ' 



R^n{l/Ts) = = -TTT^- (14) 



In the discrete-time implementation of the white noise, the sampling time At <ti Tg and Tq. 
The incoherent statistical fluctuations in the input s(i) +^(i), which controls the continuous 
noise background in the power spectral density, are measured by the variance (t| = D^/At 
|31l I33j . Here, is the rms amplitude of input noise ^{t). 

Thus, the array SNR gain, viz. the ratio of the output SNR of array to the input SNR 
for the coherent component at frequency l/Tg, follows as 

^ ^ RUl/Ts) (var[y(t)])H(l/T,)^V4' ^ ' 

Equations (|12|I - H15() can at best provide a generic theory of evaluating SNR of dynamical 
systems i31j . If the array SNR gain exceeds unity, the interactions of dynamic array of 
bistable oscillators and controllable array noise provide a specific potentiality for array signal 
processing. This possibility will be established in the next sections. 



3 Numerical results of array SR and SNR gain 

We have carried out the simulation of parallel arrays of Eq. and evaluated the array SNR 
gain of Eq. H15|l . as shown in Appendix A, based on the theoretical derivations contained in 
|31l I33j . Here, we mainly present numerical result as follows. 
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Figure 3: Numerical input SNR i?in (*), output SNR _Rout (+) at the left axes and array SNR 
gain G{l/Ts) (o) at the right function of aj, in a single bistable oscillator for (a) 

A = 0.88Ac « 0.34, (b) A = 2.QA^ « 1.0, (c) A^^iA^^ 1.925 and (d) A = 10^4^ ~ 3.849. 
Here = 100, AiAB = lO^^^ At ^ x 10~^ and = 10^ 



3.1 Improvement of the array SNR gain by noise for array size = 1 

If the array size = 1 and the response y{t) — {l/N) J2iLi Xi{t) — xi{t), this is the case of a 
single bistable oscillator displaying the conventional SR or residual SR phenomena [Tl l20[OT] . 
In Figs, m (a)-(c), we show the evolutions of E[y{t)], i?yy(T) and Cyy(T), respectively. The 
input is a sinusoidal signal with amplitude A = 0.4 and frequency l/Tg = 0.01 mixed to 
the noise ^{t). As the rms amplitude increases, the periodic output mean E[y{t)] has a 
same frequency l/Tg = 0.01, as shown in Fig. El (a), and the largest amplitude of E[y{t)] 
appears at the resonance region around — 1.1. Plots of the stationary autocovariance 
function Cyy(r), as depicted in Fig. |2 (c), indicate that the correlation time decreases 
as CTj increases, but the stationary variance Cyy(O) = (var[y(t)]) presents a non-monotonic 
behavior. As cr^ increases from 0.6 to 1.1, 1.8 and 3.0, the correlation time decreases from 
65.1 to 27.7, 17.4 and 7.5, whereas Cyy(O) equals to 0.383, 0.277, 0.390 and 0.741, respectively. 
Thus, these nonlinear characteristics oi E[y{t)] and Cyy(T) lead to the SR phenomenon of the 
output SNR i?out versus the rms amplitude crj in a single bistable oscillator, as illustrated in 
Fig. El (d) . The numerical input SNR i?in is also plotted in Fig. (2 (d) and agrees well with 
the theoretical one obtained by Eq. H14|l . Note that this SR effect is residual SR introduced 
in Ref. [201 1 since the amplitude A = 0.4 > Ac is slightly suprathreshold. Similar results are 
presented in Fig. for subthreshold amplitude {A = 0.34 < Ac) and strong suprathreshold 
ones {A = 2.6Ac, 5Ac and lOAc). Clearly, the SR-type behaviors of i?out disappear for strong 
suprathreshold amplitudes. These numerical results show recurrence of the phenomena of 
conventional SR P and residual SR [201 , and show the validity of cyclostationary analysis 
presented in SecEljSI]- 

The SNR gain G(l/Ts) is also depicted in Fig. (21 at the right axes. It is well known that 
the SNR gain G(l/Ts) is below unity, so far as the sinusoidal amplitude A is subthreshold 
or slightly suprathreshold T, 23 241 1^. as seen in Fig. El(d) and Fig. (31 (a). However, 
as the amplitude A increases to a more suprathreshold value such a,s A = 2.6Ac, G{l/Ts) 
approaches unity very closely at = 3.2, as shown in Fig. (b). Interestingly enough, 
the possibility of G{l/Ts) exceeding unity exists for strong suprathreshold sinusoidal inputs 
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Figure 4: Numerical array SNR gain as a function of the rms amplitude of array noise 
r]i{t) for (a) A = 0.34, (b) A = 0.38, (c) A = 0.4 and (d) A = 1.0. The array SNR gain curves 
correspond to TV = 1, 2, 3, 5, 10, 30, 60, 120, cxd (from the bottom up). The input noise rms 
amplitudes — 1.8 for all amplitudes A, with given input SNRs i?in = 8.92, 11.14, 12.35, 
and 77.16, respectively. Here = 100, AtAB = 10"^ At = x IQ-^ and K ^ 10^ 



{A = 5Ac or lOAc) at certain noise level regimes, as plotted in Figs. |21 (c) and (d). This 
result is consistent with the work by Hanggi et al 2 . 

These numerical results indicate that this theoretical framework of cyclostationary signal 
processing in 31 can fully describe the SR phenomena in a single bistable oscillator, and we 
shall now apply it to the SR effects in parallel uncoupled arrays of bistable oscillators with 
a noisy sinusoidal input s(t) + S,(t). 

3.2 Improvement by noise of the array SNR gain for array size N > 2 

If the array size N > 2, this is the case of parallel arrays of bistable oscillators displaying array 
SR phenomena. Figure 01 displays evolutions of the array SNR gain G{l/Ts) as a function 
of the rms amplitude a,, of array noise r]i{t), for both subthreshold {A — 0.34 and 0.38) and 
suprathreshold inputs {A = 0.4 and 1.0). The input noise rms amplitude is cr^ = 1.8, resulting 
in the given input SNRs Rin — 8.92, 11.14, 12.35, and 77.16, respectively. Then, due to array 
noise rji{t), the array SNR gain G{l/Ts) exhibits nonmonotonic behavior as a function of an 
for N > 2. This collective phenomenon can be termed as "array SR" [HSIj appearing for not 
only suprathreshold inputs, as shown in Figs. 0] (c) and (d), but also subthreshold signals, 
as presented in Figs. ^ (a) and (b). More importantly. Fig. ^reveals that the region of the 
array SNR gain G{l/Ts) raising above unity, via increasing ct,,, is possible for moderately 
large array size N. Furthermore, as A increases, G{l/Ts) reaches a larger and larger local 
maximal value for the same A^. For instance, G{\/Ts) is about 1.1 for A = 0.34 and N = 120, 
as shown in Fig. 0](a), whereas G{l/Ts) is around 1.3 for A — 1.0 and N = 120, as seen in 
Fig.il(d). 

The mechanism of conventional SR, as shown in Figs. |2HS1 exploits a combination of the 
positive role of input noise ^(f) and the nonlinearity of a single oscillator Given a 

noisy signal, the mechanism of array SR and the possibility of array SNR gains above unity 
are clearly attributed to the added array noise rii{t) interacting with the nonlinearity of the 
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Figure 5: Plots of nonstationary mean E[y{t)] of y{t) for A — 0.4 at a^j — 0, 1.3 and 4.0 
(from the top down). E[y(t)] is a periodic function of time t with period = 100. Here, 
t e [0, 2Ts[ and At = x IQ-^. 



array Figure [S] shows that nonstationary means of E[y{t)] are same for = 1, 2, 
DO, at fixed cr^, since 



N 



N 



E[yit)] = E[J2^r{t)/N] = ^i?[x,(t)]/iV = E[xS)]- 



(16) 



i=l 



However, we note that the ampHtude of E[y{t)] decreases as ct,, increases, as shown in Fig.|31 
At time t, we have 



{E[y{t)y{t + T)]) 



E 



N N 

E[x^{t)x^{t + t)] (N-l)E[x^{t)Xj{t + T)] 



N 



N 



(17) 



and 



CyyM = Ryy{T)-{E[y{t)]E[y{t + T)])=Ryy{T)- 

' E[x,{t)x,{t + t)] {N -l)E[x,{t)xj{t + T)] 



N 



N 



iV2 

E[x,{t)]E[xj{t + T)] 



(18) 



for i ^ j and i, j = 1, 2, . . . , TV. Note that E[x^{t)] = £'[xj(t)]. 

Figures|ni(a) and (c) show that, at cr^ = 1.3, Ryy{T) and Cyy(T) weaken as the array size 
TV increases. On the other hand, for a fixed array size such as = 120, Figs. El (b) and (d) 
suggest that the output behaviors of Ryylr) and Cyy(r) also weaken as cr^ increases from 
to 1.3 and 4.0. Correspondingly, the stationary variance Cyy(O) = 0.39, 0.25 and 0.15, and 
the correlation time — 17.4, 12.1 and 4.2. An association of the time evolutions of E[y{t)] 
and Cyy(r) results in SR-type curves of the array SNR gain G{l/Ts) presented in Fig. 0] 



3.3 Improvement by noise of the array SNR gain for array size 

= oo 

Figure 01 shows that the array SNR gain G{l/Ts) is an increasing function of array size N. 
Thus, it is interesting to know how much the maximal value of G{\/Ts) reaches as array size 
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Figure 6: Numerical behaviors of Ryyir) and Cyy(r). (a) Ryyir) at cr^ = 1.3 for array sizes 
A'' = 1, 2, 5 and 120 (from the top down), (b) i?yy(r) with array size N = 120 as varies 
from zero to 1.3 and 4.0 (from the top down), (c) Cyy(T) at cr^ = 1.3 for array sizes iV = 1, 
2, 5 and 120 (from the top down), (d) Cyy(r) with array size N = 120 as cr^ changes from 
zero to 1.3 and 4.0 (from the top down). Here, A = 0.4, and other parameters are the same 
as in Fig. ^ 




50 , 100 150 200 20 40 60 80 100 

time T time t 



Figure 7: Numerical output behaviors of arrays for array size N = oo. (a) i?yy(r) and (b) 
Cyy(T) as functions of ct^ = 0, 1.3 and 4.0 (from the top down). Other parameters are the 
same as in Fig. 0] 



9 



N ^oo. Form Eqs. inj and (HHJ, we have 



lim Ryyir) = lim {E[y{t)y{t + r)]) 

N-^oc N—^oo 



lim /' EMt)x,{t + t)] + {N -l)E[x,{t)xj{t + T)] 



and 



AT^oo \ N 

= Rx.xA^), (19) 

lim Cyy(r) - lim i?yy(T) - lim (i?[j/(i)]i?[y(t + r)]) 

„ , , / E[El^^^{t)]E[J:Ux,{t + r)] \ 
= ^l"?^oo^-(")-J"?^oo\ / 

= (S[a;,(t).T,(i + r)]) - {E[x^{t)]E[x,{t + t)]) 

= C,,,^.(t), (20) 

for i ^ j and i,j = 1,2, . . . ,iV. Since the indices i and j are different, but arbitrary in 
Eqs. 1)19(1 and 1)20(1 . we can adopt two bistable oscillators, each embedded with independent 
noise, to evaluate the array SNR gain of a parallel array with size N = oo. The behaviors of 
i?yy(T) and Cyy(T) arc plotted in Fig.[7|as the rms amplitude a^j increases from 0, 1.3 to 4.0. 
The stationary variance Cyy(O) = 0.39, 0.25 and 0.15, and the correlation time = 17.2, 
12.1 and 4.2, respectively. Furthermore, the output SNR i?out of a parallel array of bistable 
oscillators with infinite size = oo is obtained from Eq. 1(12(1 . and same for the array SNR 
gain G{l/Ts) of Eq. (fT^ . Numerical results of G{l/Ts) are also plotted in Fig. ^ as N = oo. 

From Figs.21-[3 the mechanism of array SR and the possibility of array SNR gain above 
unity can be explained by the fact that independent array noise, on the one hand, help 
the array response to reach its mean E[y{t)], on the other hand, counteract the negative 
role of input noise and 'whiten' the output statistical fluctuations y{t). In other words, the 
stationary autocovariance function Cyy(r) has a decreasing stationary variance Cyy(O) and 
correlation time t^, as shown in Figs. |^ and 



4 Optimization of the array SNR gain of an infinite ar- 
ray 

For a given input noisy signal and a fixed array size N, there is a local maximal SNR gain, 
i.e. the maximum value of G(l/Ts) at the SR point of rms amplitude of array noise, as 
shown in Fig. 01 Clearly, this local maximal SNR gain increases as array size N increases, 
and arrives at its global maximum G'„iax(l/T's) as = oo. Note that G'inax(l/Ts) is obtained 
only via adding array noise rii{t). It is interesting to know if G'niax(l/?s) can be improved 
further by tuning both array noise r}i{t) and the array parameter Xb- 

In Eq. 0, the signal amplitude A/Xf, ^ A is dimensionless, and the discrete implemen- 
tation of noise results in the dimensionless rms amplitude of a^/Xf, or ari/Xh Un 
(where each arrow points to a dimensionless variable). The dimensionless ratio of v4/cr^, as 
AtAB = 10~^, determines the input SNR Rm of Eq. In Fig. |H1 we adopt two given 
input SNRs i?in = 40 and 10, this is, A/a^^ ~ 0.4 and 0.2. When the array parameter X^ 
varies, but A/a^ keeps, line Li comes into being, and is divided into subthreshold region 
{A < 2/V27) and suprathreshold regime {A > 2/^27) by Hue L2 oi A = 2/V27, as shown 
in Figs.|Hl(a) and (c). We select different points on line Li, being located in subthreshold 
region or suprathreshold region, for computing Ginax(l/7s) via increasing cr^, as illustrated 
in Figs. IHl (b) and (d). Figure El (b) shows that, at the given input SNR i?in = 40, the 
global maximum SNR gain Gmax(l/7s) increases from low amplitude A — 0.25, i.e. point Pi, 
reaches its maximum around cr^ = 2.0 for A — 0.38, i.e. point P4, then gradually decreases 
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Figure 8: (a) Plots of amplitude A versus input noise rms amplitude (dimensionless 
variables). Line Li is A/a^ = 0.4, and the corresponding input SNR = 40. Line L2 of 
A = 2/\/27 ~ 0.385 divides line Li into subthreshold region (below L2) and suprathreshold 
section (over £2). Points Pi {i = 1, 2, • • • , 6) correspond to A — 0.8, 0.6, 0.4, 0.38, 0.3 and 
0.25, respectively, (b) The global maximum SNR gain G'max(l/T's), at fixed input SNR 
i?in = 40, as a function of rms amplitude ct^ of array noise for points Pi (different amplitudes 
A), (c) Plots of A versus a^. Line Li is A/a^ = 0.2, and i?in = 10. Line L2 is A = 2/^27- 
Points Qi {i = 1,2,- ••,5) correspond to A — 0.15, 0.2, 0.3, 0.4 and 0.6, respectively, (d) 
Gmax(l/rs), at i?in = 10, as a function of a,^ for points Qi. Here, = 100, AB = l/Tg and 



as the amplitude A increases to 0.8 (point Pi). The same effect occurs for the given input 
SNR i?in = 10, as shown in Fig.|Sl(d), and A = 0.2 (point Q2) corresponds to the maximum 
G'niax(l/'7s) around a,, = 1.5. These results indicate that, for a given input SNR, we can tune 
the array parameter Xfj to an optimal value, corresponding to an optimized global maximum 
SNR gain. 

However, we do not consider the other array parameter Tq, which is associated with the 
time scale of temporal variables |21| . Then, the location of optimal array parameters Xf, in 
subthreshold or suprathreshold regions, associated with optimal ct^, is pending. Immediately, 
an open problem, optimizing the global maximum SNR gain G'i„ax(l/7s) via tuning array 
parameters {Xi, and Tq) and adding array noise (increasing ct,,), is very interesting but time- 
consuming. This paper mainly focuses on the demonstration of a situation of array signal 
processing where the parallel array of dynamical systems can achieve a maximum SNR gain 
above unity via the addition of array noise. Thus, the optimization of the maximum SNR 
gain of infinite array is touched upon, and this interesting open problem will be considered 
in future studies. 



5 Conclusions 

In the present work we concentrated on the SNR gain in parallel uncoupled array of bistable 
oscillators. For a mixture of sinusoidal signal and Gaussian white noise, we observe that the 
array SNR gain does exceed unity for both subthreshold and suprathreshold signals via the 
addition of mutually independent array noise. This frequently confronted case of a given 
noisy input and controllable fact of array noise make the above observation interesting in 
array signal processing. 
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We also observe that, in the configuration of the present paraUel array, the array SNR 
gain displays a SR-type behavior for array size larger than one, and increases as the array 
size rises for a fixed input SNR. This SR-type effect of the array SNR gain, i.e. array SR, 
is distinct from other SR phenomena, in the view of occurring for both subthreshold and 
suprathreshold signals via the addition of array noise. The mechanism of array SR and the 
possibility of array SNR gain above unity were schematically shown by the nonstationary 
mean and the stationary autocovariance function of array collective responses. 

Since the global maximum SNR gain is always achieved by an infinite parallel array at non- 
zero added array noise levels, we propose a theoretical approximation of an infinite parallel 
array as an array of two bistable oscillators, in view of the functional limit of the autocovari- 
ance function. Combined with controllable array noise, this nonlinear collective characteristic 
of parallel dynamical arrays provides an efficient strategy for processing periodic signals. 

We argue that, for a given input SNR, tuning one array parameter can optimize the global 
maximum SNR gain at an optimal array noise intensity. However, another array parameter, 
associated with the time scale of temporal variables is not involved. An open problem, 
optimizing the global maximum SNR gain via tuning two array parameters and array noise, 
is interesting and remains open for future research. 
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A Numerical method of computing power spectra of the 
collective response of arrays 

The corresponding measured power spectra of the collective response y{t) ~ {l/N) J2iLi 
are computed in a numerical iterated process in the following way that is based on the 
theoretical derivations contained in The total evolution time of Eq. is {K + l)Ts, 

while the first period of data is discarded to skip the start-up transient 29, 2 . In each 
period Tg, the time scale is discretized with a sampling time At ^ Tg such that Tg — LAt. 
The white noise is with a correlation duration much smaller than and At. We choose a 
frequency bin AB = l/T^, and we shall stick to AtAB = 10"^^ ^ ^00, L = 1000 and 
K > 10^ for the rest of the paper. In succession, we follow: 

(a) The estimation of the mean E[y{jAt)] is obtained over one period [0,rs[, and the precise 
time jAt oi E[y{jAt)] (j = 0, 1, • • • , L—1) shall be tracked correctly in each periodic evolution 
of Eq. JU, i.e. [kTs, (fc + 1)T,[ for fc = 1, 2, • • • , if. 

(b) For a fixed time of r = iAt (j = 0, 1, • • • ,Tniax/Ai), the products y{jAt)y{jAt + iAt) are 
calculated for j = 1, 2, • • • , KTg/ At. The estimation of the expectation E[y{jAt)y{jAt+iAt)] 
is then performed. From Eq. ((Sj), the stationary autocorrelation function Ryy{T) can be 
estimated over a time domain r e [0, Tmax[- Immediately, the stationary autocovariance 
function Cyy{iAt) of Eq. JHl at i = 0, 1, • • • , Tt^^-k/ At can be deduced. Note the time Tmax 
is selected in such a way that at T,„ax, the stationary autocovariance function Cyy{iAt) in 
Eq. (jSJ has returned to zero. In practice, we can select a quite small positive real number e, 
such as e = 10~^. If Cyy(iAt)/Cyy(0) < e, the above computation shall be ceased and the 
index iond is found, leading to r^ax — icndAt. 

(c) Increase the total evolution time of Eq. as {K' + l)Ts {K' > K), and evaluate the 
mean E'[y{jAt)] and the stationary autocovariance function C'yy{iAt) again. If the differences 
between E'[y{jAt)] and E[y{jAt)], Cyy[iAt] and Cyy{iAt), converged within an allowable 
tolerance, we go to the next step (d). If they do not converge, the total evolution time of 
Eq. (Q) should be increased to {K" + l)Ts larger than {K' + l)Ts, until the convergence is 
realized. 
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(d) With the converged mean E[y{jAt)] and stationary autocovariance function Cyy{iAt), the 
corresponding Fourier coeflicient Yi and the power var[y(t)]_ff(l/Ts)A_B of Eq. (jHJ contained 
in the noise background around l/T^ can be nuniericaUy developed. The ratio of above 
numerical values leads to the array SNR i?out- The correlation time = AIAt as \ h{MAt)\ = 
|Cyy(MA<)/Cyy(0)| < 0.05. Thc uumerlcal input SNR ii!in can be also calculated by following 
steps (a)-(d), and compared with the theoretical value of i?i„ of Eq. The SNR gain 

G(l/T,) wiU be finally figured out by Eq. ((T^. 
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